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Abstract 

Geodesies deviation equation (GDE) is introduced. In "adiabatic" 
approximation exact solution of the GDE is found. Perturbation the- 
ory in general case is formulated. Geometrical criterion of local insta- 
bility which may lead to chaos is formulated. 



1 Introduction 



Even though classical mechanics is an old subject, the surprising fact is that 
the mechanisms affecting its dynamical evolution have only recently been 
qualitatively understood. It is known that chaos in classical mechanical sys- 
tems emerges when overlapping of resonance zones occurs. The essential 
problem in study of chaos lies in formulation of a simple mathematical cri- 
terion of sensitive dependence on initial conditions. 

An example of analytical criterion is proposed by Toda, Duff and Brumer 
[H] local criterion of the transition from the ordered motion to the chaotic one. 
The main idea of their approach lies in replacement studies of the behaviour 
of trajectories in phase space which remain nearby to the selected trajectory 
by studies of the behaviour of nearby trajectories in vicinity of selected points 
of the phase space. 

The other local criterion of chaos is studied in this article. It is based on 
observation that Euler - Lagrange equations of motion can be rewritten as 
the geodesic equation with respect to the Jacobi metric and then on mea- 
suring the local tendency of geodesies to converge or to diverge. The idea of 
this method has been originated in Krylov paper ||. The use of the devi- 
ation equation presupposes that geodesies form a congruence which is true 
locally in admissible for trajectories region of the configuration space. The 
property of being a K-system is a global property and therefore deviation 
equation which is purely local criterion needs some extension so as to obtain 
a more complete (a kind of Lapunov exponent and study of the influence of 
the boundary 8T>e on global behaviour of the system) characteristic of the 
system. 

The overall plan of this article is as follows. In section 2, we introduce 
in a new, simplified way the concept of the geodesies deviation and find a 
simple form of geodesies deviation equation. Section 3 is devoted entirely to 
searching for solutions of the GDE and to formulation of the perturbation 
expansion. The last section (4) contains examples and comments. 

2 Introduction to geodesies deviation equa- 
tion 
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2.1 Dynamics of the system in arbitrary coordinates. 

Let us consider classical mechanical system with N degrees of freedom = 
1,2..., N) described by the lagrangian 

C(q,q) = \g v q l q' -V(q), (1) 

where gijdq 1 ® dqi is a Riemannian metric on the configuration space M. 
of the system. Let us assume that the dynamics of this system could be 
equivalently described in terms of a momentum - position phase space by 
the Hamiltonian 

n{q,p) = \g ii {q)p iPi + V{q), (2) 

Pi — 9ijq denotes momentum conjugate with i—th coordinate q\ and "' = Jj" 
differentiation with respect to time. Dynamics of the system is governed by 
the Hamilton equations of motion: 

1 dg jk dV J lk 

Pi = ~Yd^ VjVk ~W q = 9 Pk ' 

By eliminating pi variables and using identity ^r9jk = ~9 %k ^f we can 
rewrite a system of the first order equations as a system of the second order 
equations 

if + r kl q k q l = -g tk d k V, (4) 
dk = and T\i is Levi-Civita connection for metric gij. 

2.2 Dynamics on the boundary of an admissible region 

An important characteristic of the mechanical system is the shape of region 
of the configuration space which is admissible for the physical trajectories. 
It is known that even small deformation of this region could lead from de- 
terministic to chaotic behavior of the system. Well known example of such 
behavior is the stadium ||. We will pay some attention to properties of this 
region. 

It is obvious that the total energy of the system is an integral of mo- 
tion. Therefore for fixed (by the choice of the initial conditions) energy of 
the system any trajectory in phase space is confined to the hypersurface 
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IgV \q)piPj + V(q) = E. The kinetic energy of the system \g l ^PiPj is positive 
and therefore (for fixed total energy E) a projection of the constant energy 
hypersurface on configuration space provides an admissible for trajectories 
(movement) region of the configuration space T>e = {? 6 M : V(o) — E}. 
T>e depending on the energy, could be bounded or unbounded, connected or 
not. In general T>e has boundary 8T>e = {q G M. : V(q) = E}. If potential 
V(q) has no critical points on the boundary then &De is Af — 1 dimensional 
submanifold of M. 

We can easy see that point q G <9£>e in which trajectory reaches 
the boundary is an isolating point in dV E . 

If trajectory reaches q G &De then speed v l in this point is equal to zero. 
It is consequence of equality \v % v 3 ' + V(q ) = E = V(q ). We see that in 
&De exist neighbourhood of qo G OT>e which does not contain any point of 
the trajectory. It means that movement along boundary 8T>e is impossible, 
and that the point go is isolating point. 

In fact, it is easy to show that if potential is smooth function of 
space coordinates then trajectories approach to q G dT>E or depart 
from q G dV E perpendicularly to the boundary dV E - 

i) Any vector £ l tangent to &D E satisfy equation 9jV(go)£ l — which 
could be rewritten in the form gij(gradV) 1 \ q = qo ^ = 0. This is just or- 
thogonality condition in the sense of scalar product defined by the matrix 
9ij(Qo)- 

ii) On the other hand if we assume that the trajectory reaches the point 
q in instant of time t then for t close to t we have q l (t) = q\ + q l (t )(t — 
t ) + ^q l (t )(t — t ) 2 + O((t — t ) 2 ). Reminding that v % = q l (t ) = and using 
equations of motion (4) we obtain q l {t) = q l — ^g i: >(q )djV(qo)(t— t ) 2 + O((t — 
to) 2 ). After differentiation of the last formula with respect to t we evaluate 
velocity in neighbourhood of the point qo, q l {t) = (gradV) 1 \ q=qo (to — t). 

From i) and ii) it follows that trajectories are orthogonal to the boundary 
dV E . 

The stadium (which is not smooth) does not satisfy assumption of the 
above statement. 
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2.3 Geometric form of the equations of motion 

Let us turn to dynamics of the system in the interior of the admissible region 
IntV E = {qeM: V(q) < E}. 

Crucial to our future construction is an observation that geodesic equa- 
tion in the Riemannian geometry defined by the Jacobi metric 
g'j = ^g 1 -" (with ip(q) = 2{E — V) and natural parameter s so that ^ = 
2(E — V)) on the set Int£> E is equivalent to the equations of motion 
(4). 

Geodesic equation for metric = ip(q)gij has the well known form 

*£ + f5*^ = , (5) 

ds 2 Jk ds ds w 

where s is a natural parameter in the sense of the Jacobi metric ( y 9ij(l( s ))%'^' = l) 
and t l j k is a Christoffel symbol with respect to the same metric If we 
express Christoffel symbols of the Jacobi metric by Christoffel symbols T l j k 
of the metric gij 

t) k = r jk + l -[d 3 {lni,)5\ + Wnifi&j - d r (lmP)g ri g jk ], (6) 

and exchange the natural parameter s for time t then we will obtain the 
rearranged equation (5) 



<*V p dq j dq k 
dt 2 jk dt dt 



ds\ 2 d 2 t dq % 
It) ds^ + ~dT 

(7) 

Now we would like to chose ip and parameter s — s(t) so as to geodesic equa- 
tion (5) be equivalent to equation (4). Let us impose the simplest conditions 
which eliminate unwanted terms from equation (7) 

( ds\ 2 d 2 t dq l n , , 



_ + —djlni, 
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Making further manipulations on parameter derivatives = %.-j t (^-jf 
we can simplify the first constrain (8) to the form 



du dlnib . . 

— = u — - — , (10) 
dt dt y ' 



where uit) = (^f) 2 - 

The solution of this equation (jj^ = u(t) = Atp 2 married with condition 
(9) gives \ diift = —diV, which provides the final answer 

ip = 2{M -V). 

M is an integration constant (A equal one was taken for simplicity). 

The only unknown is interpretation of an arbitrary constant. Reminding 
that s is a natural parameter for Jacobi Metric 

fdt\ 2 dq i dq j _ 1 d^_d^_ 
~ 9ij ds ds ~ [ )9ij [ds) dt dt - 2(M - V) 9ij dt dt ' 

we find that M is the total energy of the system 

1 dq l dq j T T ^ 
M = -ga-=--=- + V = E. 

2 y%1 dt dt 

Therefore 

tj) = 2{E-V). (11) 

We can also check that s is a good evolution parameter in the re- 
gion IntX>E i-e. t — > s(t) is monotonically increasing function of time. 

In admissible for trajectories region of configuration space %(t) = 2(E — 
V) > 0, ||(t) = -2$diV, and ||(t) = -2i?diV - 2q i qW i d j V. Let us fix 
the origin of the parameter transformation sit = 0) = 0. If for some value of 
the parameters s(t ) = s the system reaches the boundary of the admissible 
region q G dV E then 

d t fP s ffi s 

1t {to) = ~o¥ {to) = °' ~d¥ = 2 9 i3 (qo)d t V(q )d J V(q ) > 0. 

In the last inequality we used equations of motion (4). We see that so = s(to) 
is an inflection point of the considered function t — > s(t). In the region IntT>E 
s(t) is monotonically increasing function of time. 
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2.4 Origin of the GDE 



In the previous section we discovered that Euler - Lagrange equations (4) 
could be transformed to the geodesic equation 

V u u = 0, (12) 

where u is the tangent vector to the geodesic and V is the covariant derivative 
with respect to the Jacobi metric. Our purpose is to investigate the relative 
motion of geodesies in a given domain IntT>E of the configuration space. 
We choose a transverse curve C with respect to the congruence of geodesies, 
i.e., the curve which only once crosses each of the geodesies belonging to 
the congruence. We assume that the point at which the curve C crosses 
a given geodesic is the zero point of the parameter s along this geodesic. 
Now, we find the curve C s which is the copy of the curve C. We construct 
C s by transporting each point of the curve C, along the geodesic to the 
point on this geodesic at which the value of the parameter is s. In this 
way, in the domain in which the congruence of geodesies is determined, we 
obtain a 2-dimentional surface. The tangent vector fields to this surface at 
the point for which the values of parameters are (A, s) will be denoted by 
u(X,s) = |(a,s), and £(A, s) = |(a,s)- The field u is also tangent to the 
geodesic. By construction of the surface the Lie bracket of the vector fields 
vanish [u,£] = 0. Since the connection F is torsionfree, one has 

f(u,Z) = V u Z-%u-[u,£\ = 0, (13) 

and therefore 

V u £ = V^u. (14) 

Now we would like to find acceleration of the changes of the vector field £ 
along a given geodesic. From the definition of the curvature one could obtain 

R(u, f )u = (V U V ? - V ? V U - V [U)€ ]) u = V u V^u. (15) 

With the help of equations (14-15) we can find acceleration of the £ field 

V U V U £ = R(u,£)u. (16) 

The vector field £ informs us about the relative position of points on neigh- 
bouring geodesies which start from points on arbitrarily chosen curve C. We 
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are interested only in the relative motion with respect to geodesies and not in 
the motion of points along these geodesies. The component of £ parallel to the 
vector u contains useless information. Therefore, we will eliminate this com- 
ponent focusing only on orthogonal to u component of £. Let us decompose 
£ = n+Xu, where A G IR and g(n, u) = 0. We also assume unit normalization 
of the tangent field g(u, u) — 1. Under above assumption the parameter A is 
just a projection u on £, g(^,u) = g(n + Xu,u) = Xg(u,u) = X. The vector 
field n is called geodesies deviation n = £ — g(£, u)u. Since £ = n + g(£, u)u 
both sides of the equation (16) could be transformed to the form which is 
free of parallel (to the vector u) component of £ 

V U V U £ = V u {V u n + (?(V U £, u)uj = V u (V u n + g{V^u, u)uj = V u V u n, 

(17) 

R(u, ^)u— R(u, n)u + <?(£, u)R(u, u)u = R(u, n)u. (18) 

As a result of (17-18) equation (16) can be rewritten in the form known as a 
geodesies deviation equation 

V u V u n = R(u,n)u. (19) 

Equation (19) answers the question with which acceleration field n changes 
along a given geodesic. It should be stressed that equation (19) measures 
local tendency of geodesies to converge or to diverge and it works if the 
vector n is small. It also works in a neighbourhood of geodesies which are 
parallel to each other in some interval. 

2.5 Newtonian form of the GDE 

Let us rewrite the geodesies deviation equation in the form similar to the 
Newtonian equation of motion. This rearrangement allows us to use the 
known methods of solving equations of motion. 

To this end we shall use the Riemann tensor R(A, B, C, D) = g(A, R(C, D)B) 
where A,B,C,D are vector fields on configuration space M.. The compo- 
nents of this tensor are R ijM = ^(^1,^7,^,^) = 9irR r jk i- Using this 
tensor we can rearrange the right hand side of the geodesies deviation equa- 
tion (19) 

[R(u,n)uf = ~g i m klrm (u k 5 l j u r n m +u k n l u r S m j ) = ^~[R(u, n, u, n)}. 

(20) 
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Now we introduce the new gradient operator [grad n ] 1 = g % ^-§^. and define a 
"potential" V u (n) = ~R(u,n,u,n). This leads to the equation 

D 2 n - - 

— - = V u V u n = -grad n V u (n) (21) 

which looks like Newton equation of motion. 

Newtonian nature of equation (21) is even more transparent in the Fermi 
frame (E U E 2 ,..., E N _ U E N ), defined by V u E a = and g(E a , E b ) = 5 ab where 
a,b= 1,2, ...,7V- 1, N. 

If we chose the last vector of this frame as a tangent to the geodesic 
En = u and denote a, (3 = 1, 2, N — 1 then deviation equation (21) takes 
the form 

d 2 n a d , N 

~nr = — ^— V u (n). 22 
ds 2 dn a ( J ( J 

The function V is the quadratic form with respect to the n a variables V u (n) = 
A a f3n a nP, where [A a p] is (N — 1) x (TV — 1) matrix defined by the proper 
components of the Riemann tensor in the Fermi base A a p = R a N/3N- Equa- 
tion (22) describes the system of iV — 1 "coupling oscillators" with "time"-s 
dependent frequencies. If we consider the relative motion of geodesies only 
locally (in the neighbourhood of the point q* = q{s*)) then matrix [A a p] is 
constant and values of its components are defined by the values of the Rie- 
mann tensor in the point g*. It is clear that if V is positively defined then 
motion of geodesies in the neighbourhood of the considered point is stable. 
In this case geodesies approach each other. If V is undefined or negatively de- 
fined then geodesies diverge i.e. motion is unstable and we have sensitive 
dependence on initial conditions. 

Above analysis has purely local character and does not take into account 
the influence of the boundary &De on the geodesies motion. It should be 
stressed that it does not determine existence of global chaos in the system. 
Anyhow, when the admissible region T>e has no boundary then with the use 
of (22) new criterion can be formulated which determines the behaviour of 
the system completely 0. 

2.6 GDE in terms of an external curvature 

The potential V has a nice geometric interpretation. At a given point of 
the manifold, for any non-collinear vector fields A and B one can define the 
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sectional curvature Kab in the two-direction determined by these fields in 
the following way R(A,B,A,B) = K AB [g(A, A)g(B, B) - g(A, B)g(A, B)}. 
Therefore, reminding that g(n, u) — and g(u, u) — 1 we have 

1 - 1 - 

K(n) = -R(u, n, u, n) = -K un g(n, n). (23) 

Hence, the potential is proportional to the sectional curvature K un and square 
of the deviation field. 

Let us confine to the two dimensional case where (Ei, E 2 ) is Fermi frame. 
If we chose E 2 = u and n = xE\ then V u (n) = \R{u, n, u, n) = \Rivi\x 2 ■ In 
two dimensions the sectional curvature has only one independent component 
which in Fermi frame is equal to the Gauss curvature K = i?i2i2- The 
deviation equation takes now the simple form 

^ = ~K{s)x. (24) 

It is transparent that when the Gauss curvature is positive then geodesies 
approaches each other. In opposite case (K < 0) they diverge. 



3 Solutions of the geodesies deviation equa- 
tion 

3.1 Adiabatic approximation 

Geodesies deviation equation in two dimensions can be rewritten as a system 
of first order equations 

X = A(s)X (25) 
where "■" = j-. In the case of two dimensional system the vector of vari- 



es ' 

x 



ables has two components X 
A = 



x 



, and the matrix A has the form 



1 

-k(s) o 

i) At the beginning let us consider positive curvature case K(s) > 0. 
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1 





K(s) 



Having non-singular matrix V\ - 
X, X(s) = Vi(s)Y(s) and the whole equation (25) to the form 



we can transform vector 



(26) 



During this transformation only locality of the matrix V\ needs some care. 
Explicit manipulations on matrices A and V\ lead to the equation 



Y = D 1 



('- 





K 

„ 3 











Y, 



(27) 



where I is identity matrix and D 1 = V\ AV\ — y K 



-1 

1 

We see from (27) and the form of D\ that \J K has the same dimension 



as "s" derivative of K. This observation allows us to use a ratio 



K 



as a 



dimensionless perturbation parameter. 

We assume that in the neighbourhood of the considered point of the 
admissible region Gauss curvature changes very slowly in comparison to its 

value i.e. | K |<^C K2. Next we keep only zero order term in equation (27) 



Y = Di(s)y. 



(28) 



We call this simplification adiabatic approximation. Note that D\ is still 
curvature K dependent. To make further progress we denote components of 

y 



the vector Y in the following way Y 



, then equation (28) appears to 



y 

be a system of coupled first order equations 

y(s) = -^/Iis)y(s), fa) = \/X(s)y(s). (29) 

Introducing a complex variable z(s) = y(s) + iy(s) we rearrange the system 
(29) 

z(s) = i\[kz(s), (30) 

so as to easy find solution 



z(s) = z(0)exp 



ds'\IK{s') 



(31) 



11 



Coming back to the real vector Y we have 



Y{s) 



cos (j s ds'^k(s')] -sin (^J s ds'^K(s') 
sin (J, s ds'4k{7)j cos (j^ds'Jkis 7 )^ 



Y(0). (32) 



Transforming (32) by the transformation X(s) = Vi(s)Y(s) and then initial 
values of Y(0) by the transformation Y(0) = Pr 1 (0)X(0) we can express 
the vector X(s) by its initial values 



X(s) 



cos (j s ds>JIW) -jj^ sm (k S ds'4W) 



-yfk(s)sin (j Q s dsiJk(7) 



K(0) 



COS 



f s ds'Jk(t 



Therefore adiabatic problem in K > case has the solution 

x(0) 



X(0). 
(33) 



x(s) = x(0)cos ( f S ds'Jk{s')) + -^ELsin ( f S ds'Jk{s')] . (34) 



K(0) 

We see that equation (34) describes stable relative oscillations of geodesies. 
The frequency of these oscillations changes with "time" - s. 
ii) solution of the negative curvature case K < 

Construction of the solution for negative Gauss curvature is analogous to 
the proceeded construction of the positive curvature solution. The only differ- 

1 1 

ence lies in the form of the transformation matrix Vo = I — ~ — / — ~ 

L -vi ^ I vl K 

After this transformation the deviation equation takes similar to the equation 
(26) form 

Y = [V 2 - 1 AV 2 - V 2 - 1 V2]Y, (35) 



where A 



1 
I K(s) | 
icantly from equation (27) 



Detailed structure of equation (35) differs signif- 



Y = D 2 {s) /- 



I K I 

4 1 k i§ 



-1 1 
-1 1 



y. 



(36) 
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where D 2 = V 2 ~ 1 AV 2 = 



-1 
1 



For slowly varying curvature i.e. | K |<c| K |a ? we can approximate 
equation (36) by 

Y = D 2 {s)Y. (37) 

Fortunately equation (37) describes a system of decoupled first order equa- 
tions 

y(s) = -y/\ K(s) \y(s), jj( s ) = y/\ K(s) \y(s), (38) 
which has a simple solution 



y (s) = y(0)exp 



- [ S ds'J\ K(s') 
Jo 



, y( s ) = y(o)exp 



Coming back to the X vector we have the final answer 



(39) 



cosh (j s ds'^\ K(s') 



X(s) 

K(s) \sinh ( J s ds'^/\ K(s') 
First component of this vector 



sinh (j s ds'^\ K{s' 



l ^cosh(J s ds^\K(s>) 



X(0). 



(40) 



x(s) = x(0)cosh (f*ds'\J\ K(s') 



+ 



x(0) 



K(0) 



-sinh (f*ds'\J\ K(s') 



(41) 



describes exponential divergence of the geodesies with "time"-s. 



3.2 Exact solutions of the GDE 

i) Let us come back to the geodesies deviation equation (27) in the K > 
regime 

Y(s) = [D 1 (s)+e 1 (s)B 1 (s)}Y(s), (42) 

1 j tM 

and by £\{s) = -tV we denote the pertur- 



where Bi(s) 



-\ 



bation parameter. 

To define perturbation procedure we introduce a kind of interaction pic- 
ture. Having exact solutions of the adiabatic problem Y(s) we define vector 
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U(s) which carry part of the evolution of the system which goes beyond the 
adiabatic approximation solution 



Y(s) = 



cos (j s ds'^K(s')) -sin (^ s ds'Jk(s 



sin J s ds'^JK(s') cos J s ds'^K(s>) 



U(s). (43) 



With the use of the new vector field the equation (42) can be rewritten in 
the form 

U(s) = ei(s)Ci(s)t/(s), U(0) = Y(0), (44) 

where 



Ci(s) = --y/K 



The formal solution of the problem 



sin 2 (f s ds'JI(^)) \sin (2f s ds>JI(7) 



\sm 2f s dsyK(s') cos 2 j s ds'^K{s' 



U(s) 



1 + 



rs rsi [Sk-i 

> / dsi ds 2 ... / ds k e 1 (s 1 )e 1 (s 2 )...e 1 (s k )C 1 (s 1 )C 1 (s 2 )...C 1 (s k ) 

k=l JQ J JO 



Y(0) 



(45) 

could be treated as a definition of perturbation theory. If we cut the series 
on arbitrary k = N then we obtain approximate solution of the problem 
(44). The precision of cuted solution (45) is only limited by the number of 
considered terms N. 

ii) In the negative Gauss curvature regime geodesies deviation equation 
has the form 

Y(s) = [D 2 (s) + e 2 (s)B 2 (s)]Y(s), (46) 

-1 1 
1 -1 



where B 2 (s) 



and perturbation parameter e 2 (s) 



\ K }($ . Now we transfer a part of the time dependence of the Y(s) on vector 
\k\* 

U (s) which carries non-adiabatic part of the evolution 



Y(s) 







cxp (-/o'daVl K(s') 

exp[J s ds^\K(s') 



U(s). (47) 
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Hence U(s) satisfy 



where C 2 (s) = \y/\ K 
The solution in this case 
U(s) 



U(s) = e 2 (s)C 2 (s)U(s), U(0) = Y(0), (48) 
-1 exp (2f a ds'y/\ K(s') 



exp 



-2f '<tfJ\ K{s>) 



-1 



U(s). 



I +V/ dsi \ ds 2 ... ds k e 2 (si)e 2 (s 2 )...e 2 (s k )C 2 (s 1 )C 2 (s 2 )...C 2 (s k ) 
k=1 Jo J o Jo 

(49) 

could also be used to construct the perturbation expansion. 



no), 



4 Stability in examples. 

In further discussion we will use formula which express Gauss curvature 
through total energy of the system, potential and its derivatives. Let us 
find explicit relation between the Riemann curvature tensor R\- k for Jacobi 
metric and curvature tensor R\jk for metric 



R l jki — R\m + 



jki 



N-2 



[5 l iCjk — S\Cji + gkjC l i — gjiC l k], 



(50) 



where auxiliary tensor CV, is built of total energy and potential in the fol- 
lowing way 



(N-2) 

13 A{E-vy 



2(E - V)ViVjV + 3(V i F)(V i y) - -g l3 M v 2 



N is number of dimensions and A/y 2 is square of the gradient of the potential 
jVy 2 = g tJ (diV)(djV) which is positive quantity. 

Equation (50) allows to express the curvature scalar for Jacobi metric 
through the curvature scalar for metric 



R = g ij R i3 = g ij R k ik3 



-R^ 



(N-l) 



2(E - V) 8(E - Vf 



4(E - V)AV + (6 - N)j\f\ 



v 



(51) 
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where R = g %3 Rij = g l ^R k ikj and AV denotes Laplacian AV = g^ViVjV. 

Large number of mechanical systems is determinated by metric gij which 
is just flat metric transformed to curvilinear coordinates. If we confine our 
interest to this class of systems then first term in equation (51) disappears 
[R = 0). Let us consider two dimensional case (N = 2) where exists ex- 
tremely simple relation between scalar curvature and Gauss curvature of a 
given surface K = ^R. Equation (51) in two dimensions leads to simple form 
of Gauss curvature 



K = 



4(E-vy 



(E-V)AV + M V 2 ] . (52) 



In considered region of configuration space IntT>E we have positive (E — V) > 
and the only term which can affect the overall signe of the Gauss curvature 
is the Laplacian of the potential AV. 

Now we are ready to test criterion (52) on two dimensional systems. 

Example 1 

Let us consider lagrangian which describes small oscillations 

C = i(i 2 + y 2 ) - {ax 2 + 2(3xy + iy 2 ). 

Potential in this case is quadratic form of coordinates. We can check that 
if V is positively defined (i.e. a > and > f3 2 ) then Laplacian of the 
potential is positive AV — 2(a + (3) > and therefore K > 0. According to 
our criterion (52) trajectories do not diverge. 
Example 2 

Another example is two - body problem which is defined by the lagrangian 

C = ^(r 2 + r 2 <p 2 )-V(r), 
where a = mi , m2 is reduced mass of the system consisting of two masses 

™ mi+rri2 J 

mi, m 2 separated by the distance r. This problem is integrable and angular 
momentum is its integral of motion L = fir 2 Lp. 

Jacobi metric for this system is equal = 2{E — V(r))gij, where [g^] = 

// ^ 2 . Gauss curvature of the geometry fixed by Jacobi metric is 



r 2 
K 



Afi(E - Vf 



(E - V) (V"{r) + l -V\r)) + {V\r)f 
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For particular form of the potential V(r) = —-, (a > 0) Gauss curvature is 
even more simple 

k(r) = ^ . 

In the Kepler problem we have two cases. 

i) We have eliptic motion in the regime E < . In this case K > and 
trajectories do not diverge. 

ii) If total energy is positive E > then K < and we have hiperbolic 
motion. Although trajectories diverge the system is integrable. 

The exponential instability (i.e. sensitivity to the initial conditions) is one 
of the important features of the chaotic evolution of dynamical systems and 
therefore determining of the simple and possibly elegant criterion decisive 
when this instability will take place, is so important. 

Our formula (22) and the comments to it, show the elegant, geometric 
formulation of the criterion of the local stability of relative movement of 
geodesies. The chaotic behaviour of dynamical systems is very complicated 
phenomenon and the exponential instability (the negative sign of Gauss cur- 
vature K) is not mostly the sufficient condition in order to the corresponding 
systems will behave chaotically (as the example of Kepler problem shows). 
It is known from Anosov paper || that on the compact manifold without 
boundary the geodesies family corresponding to the geometry with everywhere 
negative Gauss curvature K behaves chaotically. Unfortunately, space T> E 
admissible for movement, usually has non - empty boundary 8T>e, on which 
the corresponding geometry, determined by Jacobi metric [g~ij] becomes sin- 
gular and which has important influence on the geometry of the family of 
geodesies. Analysis of this influence requires application of more advanced 
global methods. 
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